The Approximate Invariance of the Average Number of Connections 
for the Continuum Percolation of Squares at Criticality 
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We perform Monte Carlo simulations to determine the average excluded area < A ex > of randomly 
oriented squares, randomly oriented widthless sticks and aligned squares in two dimensions. We find 
significant differences between our results for randomly oriented squares and previous analytical 
results for the same. The sources of these differences are explained. Using our results for < A ex > 
and Monte Carlo simulation results for the percolation threshold, we estimate the mean number 
of connections per object B c at the percolation threshold for squares in 2-D. We study systems of 
squares that are allowed random orientations within a specified angular interval. Our simulations 
show that the variation in B c is within f .6% when the angular interval is varied from to 7r/2. 



I. INTRODUCTION 

Continuum percolation has been of significant inter- 
est in the study of porous media [[l] . It offers important 
advantages over lattice percolation due to the fact that 
the majority of systems encountered in nature are not 
confined to a lattice and are therefore modeled more ap- 
propriately using continuum systems ■ 

When studying the transport properties in porous me- 
dia, the connectivity properties of the spanning cluster 
at percolation threshold are important. One measure 
of the connectivity is the mean number of connections 
per site. In the case of lattice percolation, Scher and 
Zallen || demonstrated the approximate dimensional in- 
variance of this quantity. The behavior of the analogous 
quantity B c in continuum percolation systems has been 
previously studied || |). In the case of continuum per- 
colation, the product of the critical concentration N c of 
objects at the percolation threshold and the average ex- 
cluded area < A ex > gives the critical average number of 
intersections per object B c 



B r = Nr. < A K 



> 



(1) 



The excluded area of an object is defined as the area 
around an object into which the center of another simi- 
lar object is not allowed to enter if intersection of the two 
objects is to be avoided |9j. In the case of objects that 
are allowed random orientations in a specified angular 
interval, one defines an average excluded area < A ex > , 
that is the excluded area averaged over all possible ori- 
entational configurations of the two objects. It has been 
claimed J5J that B c for percolating systems of differently 
shaped objects lies within a bounded range : in 2-D, 
3.57 < B c < 4.48. B c represents the connectivity in the 
spanning cluster and is of interest as the invariance of B c 
would enable us to estimate the percolation threshold N c 
using Eq. ([j]), once < A ex > has been calculated. 

In the present work we focus on continuum percola- 
tion systems of squares in 2-D, in which the objects are 



allowed random orientations within a specified angular 
interval. The motivation for the study of such orienta- 
tionally constrained systems comes from the geological 
observation that fractures in rocks do not have random 
isotropic orientations but are oriented within a more or 
less fixed angular interval. For our system of squares 
there is one angle that specifies the orientation of the 
object and we constrain it to lie within —9^ < 9 < 9^. 
We determine the percolation thresholds for different val- 
ues of the constraint angle 9^. Simulations are also per- 
formed to find the excluded area for each case. Our re- 
sults show that for a given object shape B c is constant to 
within 1.6% for squares independent of the value of the 
constraint angle 9^. 



II. SIMULATION METHOD FOR FINDING 
EXCLUDED AREA 



Here we describe the method used to determine the ex- 
cluded area for a pair of objects that are allowed random 
orientations within the angular interval —9^ to 9^. For 
rectangles (squares being a particular case) 0^ = corre- 
sponds to the case where the objects are aligned parallel 
to each other and 9^ = ir/2 corresponds to the random 
isotropic case. We describe the algorithm for finding the 
excluded area of squares in 2-D. 

A square of unit side is placed with its center coin- 
ciding with the center of a lattice of edge length L = 5. 
The lattice size is chosen to be larger than the excluded 
area, but small enough to sufficiently minimise the num- 
ber of wasted trials and yield good statistics. The square 
is given an orientation 9i, randomly chosen in the inter- 
val — 9^ < 9i < 9^, with the reference axis. A second 
square is then introduced into the lattice with its center 
randomly positioned in the lattice. This square is given 
an orientation 0j chosen randomly from the same inter- 
val as for the first object. We then determine if the two 
squares overlap. We repeat this procedure for 10 9 trials 



1 



and record the number of times the two squares overlap. 
This number divided by the number of trials is the prob- 
ability that the two objects overlap. The probability of 
overlap times the area, L 2 , of the lattice yields the aver- 
age excluded area for a pair of squares oriented randomly 
between —9^ to 9^. The method used to determine the 
overlap of squares in 2-D is described in detail in Ref. 




III. EXCLUDED AREA SIMULATION RESULTS 

We determine the average excluded area for a unit 
square for different constraint angles. We also determine 
the excluded area of widthless sticks for the case of ran- 
dom isotropic orientations. In our simulations the width- 
less stick is represented by a rectangle of edge lengths 1 
and 1 x 10~ 12 . Table 1 lists the Monte Carlo results for 

< A ex > obtained for the different objects studied for 
the case of random isotropic orientations. Table 2 shows 
the variation of < A ex > for squares with the constraint 
angle. Our values for < A ex > for aligned squares are 
consistent with all previous results 0. Furthermore, our 

< A ex > values for randomly oriented widthless sticks 
are also consistent with earlier determinations 0; our 
slightly higher value of < A ex > compared to that of Ref. 

is explained by the fact that our widthless sticks have 
a finite width and thus are expected to exhibit a larger 

< A ex > than found analytically for the zero width limit. 
However our < A ex > for the case of randomly aligned 
squares is different from previous analytical results, our 
value being 12% above that determined by Ref. 0. We 
propose a reason for this difference in the next section. 



IV. DISCREPANCY WITH PREVIOUS 
ANALYTICAL DETERMINATION OF 
EXCLUDED AREA 

We investigated the cause of the difference between our 
value of < A ex > for squares in 2-D and the previous an- 
alytical result in and found it to be the following: In 
arriving at the expression for < A ex > Ref. |?J finds the 
excluded area for a pair of rectangles (Eq. 18) with a 
given relative orientation 9 (see Figure |l|). For squares, 
using equation (Eq. 18) in Ref. |7j], 

A ex = (sin 9 + cos 9 + l) 2 - 2 sin 9 cos 9, (2) 

where 

9^\9 % -9 3 \, (3) 

9i and 6j being the individual orientations of the two 
squares. Ref. then obtains the average excluded area 
by averaging the right hand side of Eq.(||) over all pos- 
sible orientations of both objects, —9^ < 9i < 9^ and 
—Oft < 9j < 9^ , using a uniform probability distribution 



P(9 l ) = P(6j) = 1/2^. (4) 

However, it appears Ref. overlooked the fact that 
Eq.(|) holds only for < 9 < tt/2 (hence < M < tt/4), 
since for 7r/2 < 9 < tt, the expression gives a value of A ex 
less than the minimum possible value of 4 0. Thus the 
procedure of Ref. does not work when the constraint 
angle 6* M is greater than tt/4. The correct result can be 
obtained by replacing 9 in Eq.(Q) by 

9' = 9 mod (tt/2), (5) 

so that Eq.(||) holds for all values of For the ran- 
dom isotropic case 6> M = tt/2, using the modified Eq.(j|) 
and integrating numerically we obtain < A ex >= 4.54647 
which is in close agreement with our Monte-Carlo simu- 
lation result. 

We also calculate the values of < A ex > for squares 
with other values of 9^ between and tt/2 (Table 2). 
We notice that the values of < A ex > are the same for 
9^ = Tr/A and 9^ = tt/2, which is true since the rota- 
tion of a square in a particular configuration through an 
additional angle of 7r/4 yields the same configuration. 

Note that the value of < A ex > appears to decrease 
for 9^ > tt/4: , reaches a local minimum near 9^ = tt/3 
and then increases again till it reaches a maximum at 
Of! = tt/2 (see Fig. ||). This can be explained as follows. 
The case M — tt/4 is equivalent to the case of random 
isotropic orientation. Here the angle 9 — \9t — 9j\ can 
range from to tt/2. When 9^ is greater than tt/4, 9 can 
take values greater than tt/2 which means that in addi- 
tion to the configurations obtained for # M = 7r/4, there 
are other configurations for which the relative orientation 
tt/2 < 9 < 29^. However the latter configurations are, 
in fact, the same as those for 9 = (29^ — tt/2) < tt/2 
due to the symmetry of squares. Thus, the decrease in 
the < A ex > between # M = tt/4 and 6* M = tt/2 can be 
attributed to the increased probability of achieving con- 
figurations with smaller excluded areas. 

V. PERCOLATION THRESHOLDS 

Using the procedure of Ref JT3], we perform Monte 
Carlo simulations for the determination of the percola- 
tion threshold based upon the Leath method [|ll| and 
the methods Lorenz and Ziff |l2| used in their study of 
continuum percolation of spheres. The only difference 
in our present simulations is that the random numbers 
generated to fix the orientation of an object lie within a 
specified angular range from — 9^ and 9^. We determine 
the percolation thresholds of squares in 2-D for differ- 
ent values of the constraint angle 9^. These results are 
shown in Table 2. We also show the values of critical area 
fraction (f> c . Fig. || shows a plot of the percolation thresh- 
old N c for squares in 2-D versus the constraint angle 0^. 
An interesting feature of the 2-D plot is that as we in- 
crease the orientational freedom beginning from 9^ = 0, 
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N c drops until it reaches the value for 9^ = 7r/4, then 
begins increasing until it reaches a maximum and then 
falls again. This behavior is expected if B c is to remain 
approximately invariant, as we shall explain below. 



the value B c can be considered as a unique quantitative 
characteristic of a shape and can therefore be useful in 
the prediction of the percolation threshold as has been 
previously pointed out j| |j| . 



VI. APPROXIMATE INVARIANCE OF B c 
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Using the percolation thresholds and excluded area ob- 
tained from our Monte Carlo simulations, we find the av- 
erage number H of connections per object at threshold 
B c . Table 2 shows the values of B c for squares for various 
values of the constraint angle The change in B c in 
going from a constraint angle of 9^ = to # M = tt/2 is less 
than 1.6%. Using the old values of < A ex > @, the vari- 
ation in B c is seen to be rj 9.5%. We see that the slight 
decrease in < A ex > between 9^ = ir/4 and 6^ — tt/2 (see 
Fig. |^) is compensated by an increase in the correspond- 
ing iV c (see Fig. ||) to give an approximately invariant B c . 
The closeness of B c values for a given system is striking 
and is consistent with the hypothesis that B c is approxi- 
mately invariant for continuum percolating systems of a 
particular shape. 



VII. SUMMARY 

Our results show that the value of B c is approximately 
independent of orientational constraints. The B c value 
of a shape is indicative of the efficiency of the object 
in forming a percolating cluster. Not only the magni- 
tude of the excluded area plays a part in the formation 
of connections, but also the distribution of the average 
excluded area in space. This is easily seen from the fact 
that both unit area discs and aligned unit area squares 
have < A ex > = 4 j7|, but the percolation threshold 
of aligned squares <f> c — 0.6666 [|10[ is lower than that 
of discs 4> c — 0.676339 p3]. Our results suggest that 
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TABLE I. Comparison of the average excluded area for widthless sticks and squares in 2-D. The uncertainty in < A ex > 
is estimated as follows. The reciprocal of the square root of the number of Monte Carlo trials yielding intersection of the 
two objects is the fractional uncertainty in the determination of < Aex >■ The product of the fractional uncertainty and the 
estimated value of < A ex > is the uncertainty in that value. 



Object 


< Aex > for unit area object 


Previous Result 


Widthless sticks 


0.6367 ±0.0001 


0.6366 § 


Aligned squares 


3.9998 ± 0.0003 




Random Squares 


4.5466 ± 0.0004 


4.084 g 
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TABLE II. Critical area fraction, percolation threshold, average excluded area and critical average number of connections 
per object for squares in 2-D. Estimation of uncertainty in 4>c is described in E3|. 



Constraint Angle 


4>c 


N c 


< A ex > 


Be 


tt/2 


0.6254 ± 0.0002 


0.981896 


4.5466 ± 0.0004 


4.464 


tt/3 


0.6265 ± 0.0005 


0.984837 


4.5309 ± 0.0004 


4.462 


7r/4 


0.6255 ± 0.0001 


0.982163 


4.5459 ± 0.0004 


4.465 


tt/8 


0.6355 ± 0.0005 


1.009229 


4.4076 ± 0.0004 


4.448 


tt/16 


0.6485 ± 0.0005 


1.045546 


4.234 ± 0.0004 


4.443 


tt/32 


0.6575 ± 0.0005 


1.071484 


4.1240 ± 0.0004 


4.419 





0.6666 ± 0.0004 


1.098412 


3.9998 ± 0.0003 


4.394 
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FIG. 1. Procedure for determining the excluded area of two squares of side L: the first square i (shaded), is kept fixed while 
the second square j having orientation 8 with respect to i, is moved around i always keeping contact, and the locus of the center 
of j is found. The area within the locus gives the excluded area for a given relative orientation of the two squares. 
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FIG. 2. Plot of the excluded area < A ex > for squares in 2-D versus the constraint angle M . The uncertainties are smaller 
than the symbols. 
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FIG. 3. Plot of the percolation threshold N c for squares in 2-D versus the constraint angle M . The uncertainties are smaller 
than the symbols. 
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